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ABSTRACT 
Two new transient testing techniques 
method developed by Kohlmayr and the time 
zero intercept method was found to be the 
is limited to values of Neu < 2.5. The 


effectively when the value of Neu is less 


were evaluated; the centroid 

zero intercept technique. The 
most promising of the two but 
centroid technique can be used 


than 5.0. 


A heater system made of .001 inch diameter nichrome wire was designed 


and tested to determine its effect on the 


type heat exchangers. Because the design 


transient testing of matrix 


showed no improvement in the 


test results and was unreliable its use was discontinued. 
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NOMENCLATURE 


English Letter Svmbols 


A 


A 
Cc 


Matrix total heat transfer area 
Matrix iiniin- free flow area 
Matrix total frontal area 


Matrix solid cross-sectional area 
available for thermal conduction 


Fin thickness 
Flow passage perimeter (A,) 


Fluid stream thermal capacity rate 
(me ,) 


Fluid specific heat 
Matrix thermal capacity (Wie. 


Fluid specific heat at constant 
pressure 


Matrix material specific heat 


Flow passage hydraulic diameter 
(4r, ) 


Friction power per unit area 
Flow stream mass velocity (m/A_) 


Proportionality factor in Newton's 
Second Law 


Normalized fluid temperature at 
inlet to test section 


Surface heat transfer coefficient 
for convection; heat transfer 
power per unit area per degree 
temperature difference 

Deviation from step 


Fluid thermal conductivity 


Matrix thermal conductivity 


sq ft 
sq ft 
sq ft 


sq ft 


mt 
ft 


Btu/(hr deg F) 


Btu/(lbm deg F) 
Btu/deg F 


Btu/(lbm deg F) 


Btu/(lbm deg F) 


ft 


hp/sq ft 
lbm/(hr sq ft) 


32.2 (lbm ft)/(1lbf 
sec“) 


dimensionless 


Btu/(hr sq ft deg F) 


dimensionless 
Btu/(hr sq ft deg F/ft) 


Btu/(hr sq ft deg F/ft) 


NS 


Total matrix flow length 
Mass flow rate 

Pressure 

Matrix porosity (A. /A,) 
Heat transfer rate 

Gas constant (53.35-air) 
Hydraulic radius (AL/A) 
Temperature 

Normalized temperature 
Flow velocity 

Matrix volume 

Matrix mass 

Fluid mass in matrix 


Distance along flow passage 
from the matrix inlet 


Reduced length (N_ x) 
tu L 


Greek Letter Symbols 


Compactness (A/V) 


Ratio of orifice diameter to pipe 
diameter (do/d) 


Difference or change (time, 
temperature, distance, etc.) 


Time 

Free time 

Fluid viscosity 
Centroid coordinate 


Density 


10 


fe 

lbm/hr 

lbf/sq ft 
dimensionless 
Btu/hr 

(ft lbf)/(lbm deg R) 
i 

deg F 
dimensionless 
ft/sec 

Cuore 

lbm 

1lbm 


fat 


dimensionless 


sq ft/cu ft 


dimensionless 


sec, hr 
dimensionless 
lbm/hr ft 
dimensionless 


lbm/cu ft 


Subscripts 


atm 


ave 


f 


Local atmosphere 

Average 

Fluid (gas, air) 

Initial, inlet 

Matrix, mean 

Ae orbfice 

Solid (Matrix material), static 
Standard (temperature and pressure) 
Local conditions 

Inlet conditions (upstream of matrix and heaters) 
Inlet conditions at matrix entrance 


Exit conditions at matrix outlet 


Dimensionless Groupings 


f 


Fanning friction factor; ratio of wall shear 

stress to fluid dynamic head 

: 2/3 

Colburn j-factor (N Np ). This factor plotted vs. 
Reynolds Number de fities the surface heat transfer 
characteristics. 


Longitudinal heat conduction parameter for solid 
material (k A. /mLc -) 


Time parameter (hAG/W ic) 


le 
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ie, Introduction. 

The transient test facility at the Naval Postgraduate School (NPS) 
has been in operation for several years. Two techniques have been used 
at the NPS facility for the determination of heat transfer data. One of 
these is the maximum slope technique developed by Locke [13] which uti- 
lizes the temperature-time response curve of the fluid leaving the matrix 
after a step change in the fluid inlet temperature. Howard [6] extended 
this technique to include the effects of longitudinal conduction. The 
other technique is the cyclic technique developed by Bell and Katz [4]. 
In both of these techniques the heat transfer parameter being sought is 
Neu? which is a dimensionless heat transfer parameter equal to the ratio 
of the convective heat transfer rate from a solid to the heating capacity 


rate of an adjacent fluid, i.e., 


) ZA 
where: PRE 


PD = unit conductance for convection heat 
transfer (BTU/hr sq ft deg F) 


Af = total heat transfer area of solid (sq ft) 
mM = mass flow rate of fluid (lbm/hr) 
XY = specific heat of fluid (Btu/lbm deg F) 

It is known that the maximum slope technique is unreliable at values 
of Neu less than 3.5, due to the large errors in Ms associated with 
errors in the determination of maximum slope. Furthermore, it has been 
noted experimentally that for values of Neu near 2.0, the temperature 
response curve of a fluid displayed its maximum slope at approximately 
time zero on the trace, and that this value of maximum slope was not 


suited for determining Ras 
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It was not until the work done by Kohlmayr [9, 10], which gave the 
exact analytical solution to the single-blow problem based on Hausen's 
[5] mathematical model, that the reason for these errors was fully under- 
stood. Later work by Kohlmayr [11, 12] demonstrated ways in which to 
handle this problem as well as the development of a new technique which 
might be used in place of the maximum slope technique for values of Ney 
<5. Kohlmayr also developed a means to handle the single-blow problem 
when other than a step change in the inlet fluid temperature was made. 

The purpose of this thesis was first to determine experimentally 
the actual inlet temperature response of the NPS facility and to use this 
known temperature response in the manner suggested by Kohlmayr. Addi- 
tionally, a new heater system was made to try to more closely approxi- 
mate a step change in the fluid inlet temperature so that the physical 
process and mathematical model might more closely resemble one another, 
thereby improving the results using the maximum slope technique. 


Another technique, the “time zero intercept technique," used by 


Wheeler [17] was investigated. 
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2 Summary of Theory. 
A. Background 

The single-blow transient technique which is used to determine 
heat transfer data for a porous solid originally used for its mathe- 
matical model Hausen’'s [5] partial differential equation system. The 
method involved comparing the recorded exit temperature of a fluid passing 
through a porous solid which had previously undergone a step change in its 
temperature, with a computed response curve based on the solution to Hausen's 
equations. 

Locke [13] has shown that there exists a unique relationship 
between the maximum slope of the response curves and the number of heat 
transfer units, Neus However, Hausen's model did not include the effects 
of longitudinal conduction. Howard [6], by the use of a finite differ- 
ence technique with the digital computer, included the effects of longi- 
tudinal conduction. 

It should be noted here that in all of Kohlmayr's work, which 
will be discussed later, the effects of longitudinal conduction are not 
considered. 

B. Theory 

The basic assumptions in the single-blow problem are: 

(1) Properties of the fluid are temperature independent 

{2) Fluid flow is steady 

(3) The porous solid is homogeneous 

(4) The thermal conductivity of both solid and fluid is 
infinite in the direction perpendicular to flow 

(5) Thermal conductivity of the solid is zero in the 


direction of flow 
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Hausen's original differential equation system [5] was based on 
the energy balance between a fluid passing through a porous solid and the 
solid. The system of differential equations which resulted from this 
development are as follows: 

It Cr ens 
socllaiiiec ze = 
oO x 
At es 
aes 
At 
in which a "reduced length" variable, 
z= thx -N% ox 
Li goulk “oe 
and a dimensionless time parameter, 
C= BA a AA Wye A¥ 
We, & 1) L. Le \W 4s 


have been used. 


oe 





In the above equations: 


ec temperature of the solid (deg F) 
ty = temperature of the fluid (deg F) 
© = time (hr) 
Xs = specific heat of the solid (Btu/lbm deg F) 
Wis = mass of the solid (1bm) 
SX = distance along flow passage measured from inlet (ft) 
lL = total length of solid (ft) 
Wy = mass of fluid entrained in solid (l1bm) 
Kohlmayr [11] has modified Hausen's original equations by introducing a 


new dimensionless time variable called "free time," 


LF = NLF a= We ke (1) 
s Xs Ws Ls 


or 
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_". 
a, e "= 


For convenience Kohlmayr introduced two constants, 


ea = TOA and f= Wye L&E 
Ws Ls sks 


which are fixed for any particular experiment. For most practical pur- 


poses, G0 and 4 | 0s 


Hausen's modified equations are then given as 
dtlen) le”) _ ere 
GZ 


It (ZA) EM) ma Po), 
Neu Of 


(3) 
(4) 


whe re tz and te are the temperatures of the solid and fluid respec- 
tively, which have arbitrarily been normalized. The above equations are 
subject to the boundary and initial conditions: 


at Z = 0 Ce Bw) = HY) (5) 


and 
at 4 = 0 BPEewM) = | (6) 


where g( X ) represents the normalized, time-dependent fluid tempera- 
ture at the inlet cross-section. Kohlmayr solved equations (3) through 
(6) by means of a double Laplace Transform [10]. 


The results of —— s solution are: 


tole w) = (8M) [2M = 2M, Gen) EK “andv) 


Cis 
+1-e [i+ J zN..=. (2N,,V) eMYy vy] 


—\ is an entire family of functions of the order K (K=0,1,..,1N) 


i.e = 
cr ee 
nNiMakyi 
Evaluated at N=O 
a ie > Ne (8) 


equation (7) becomes 


iil 
OF Mey ft) = I EC * |) GY? 


(9) 
ye 
Se ~Me, (Vf) 
i IX: =, (M2 vw)e Me VA? 6, _ WVIAV 
By differentiation of equation (9) with respect to 44 ; 
“MulIFUM » 2 MM em Une 

dt (M,. ) = e oe, GUA) — Nee ar Li, JA.) 

on i. 
(10) 


ph 
53 — te Zz -NM,,, (U-V) 
a Ime A =O (u-v)) Te, (WM, “ue-n) | eo wu Wy) AV 


With these results Kohlmayr [12] then showed that for other 
than a step change in inlet temperature, both the maximum slope and "free 
time" at a given value of Neu may be multi-valued. He further demonstrated 
that by knowing what the inlet response is, one can determine which of 
the values of maximum slope is valid and also how to determine what the 
error bounds are in the use of some known inlet temperature change. One 
of Kohlmayr‘s results is that the correct value of "relative" maximum 


slope occurring in the fluid exit temperature response will be that value 
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which occurred last in cases in which Kohlmayr's technique is applicable. 
Furthermore, the maximum slope method is unstable for Neu < . 2.0. de 
singular at Neu = 2.0, and inaccurate for Pea between 2.0 and 3.0. [9,12] 

In view of the short-comings of the maximum slope technique, 
Kohlmayr developed an indirect curve-matching technique [11]. Some im- 
portant results of this development are summarized below. 

With the solution to the single-blow problem known, it is pos- 
sible to develop an indirect curve-matching technique based on the first 
moments of the fluid transient response curves. This method, known as 
the centroid method, involves the reducing of both the theoretical tempera- 
ture response te(N i 4) and experimental response ee into 
two different single-valued functions based on the one parameter Neu! 

To generalize the problem further, a mapping functional was defined, 
based on the fluid inlet temperature change and the fluid's exit tempera- 


ture response such that 
— 
A (Nu } a srw, wee 9 QU) 


where the following restrictions are imposed on V: 
(1) VY“ must be real, single-valued, continuous with 
> ~— 
respect to both ce(N, ek ) and g( AL), and monotone with the parameter 


N e 
cu 


(2) y 5 i must be monotone increasing with I(g) = 


[aude , the "deviation from step.” 
a LN) 
=P CLNeu 


for the maximum permissible amplification of errors. 


> 
a , where K is some measure 








(4) For any Neu and any given deviation from step,  (N,) 


should be insensitive with respect to local variations of g( Ac). 


(5) The evaluation of Ve (ND must be simple and 
straightforward. 
The functional chosen was the first moment of the difference 


between upstream and downstream fluid temperatures, 


[He [6 has) ~ ashe = b+ b+ [alge ay 


Due to the difficulty involved with integrating equation (11) up to 


values of {4 =@ , anew functional was chosen and defined as 


Jp [Ea -942)| be = Za + [Hence (12) 


Some further restrictions imposed upon g( LL) are; 
(1) eC LL ) must be non-negative monotone decreasing 
(2) have initial value g( O ) = 1 and 


(3) assume value zero for all free times which exceed 
the maximum permissible deviation from step: gC ir ) 


= 7 
0 for pict cs Tenge 0 


Therefore the moment functional was defined in terms of the centroid of 


the area under this difference curve: 


L y ti (ALP) - pe) | hee a 


LEENtb i. 14, 
— Ma) — ge) | be 
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For the case of eC LL ) = 0 (step change), the results of equation (13) 
are listed as Figure l. 

In order to use the above results for any given inlet tempera- 
ture change BC AL ), Kohlmayr arrived at the following empirical relation- 


ship: 
hese Men 9) Hsp ( Mn 19) + 99 Ly) 4 (14) 


This relationship was obtained by systematically applying the centroid 
method to many different upstream fluid temperature changes. 
Kohlmayr's solution to the single-blow problem may also be used 


directly. Returning again to equations (9) and (1) 


—N 


67 (Mw) = I- ee ™* | 1-9) 


jt . rls, (Vf) 
(M22, OL) € C9 Vv) LY 


where 
HE =—ak Me ie 7 ) Wy Xe Tk 
W- ¢, | es 
tx(N ) can be evaluated at "free time," < = 0, which would corres- 
pond to the time that an element of fluid which has undergone a change in 


temperature upstream arrives at the exit of the solid. Thus, 


=H 
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) ~ tio) = &%& [1-900] 


(16) 


Taking the natural logarithm of both sides of the above equation, yields 
— ” —~— aoe ; — 
Ly / as iW, ,0)| = M+ 4n| | go) it 
Therefore 
fee _ 4* ig 
Noa Le [ tfcy,,,0)| + 4n / Glo) (18) 
Note that for a step change in the inlet temperature, i.e., g( KX) = 0 
for A >0, that equation (18) reduces to, 
= 8 sg Se 
N= Wr [/ ty NO (19) 


In order to utilize equation (19) directly it would be necessary to achieve 
a step change in the inlet fluid temperature, g(r ). This would occur 
at real time, © =0. Then measure the temperature of that same element 
of fluid which had undergone the step change in temperature, as it left 


the porous solid. This would occur at "free Sse = 0. 
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Because of the finite response time of the different components 
of the experimental equipment, i.e., heaters, thermocouples and recorder, 
it is not possible to meet precisely either of the above requirements. 

In view of the above limitations it is necessary either to make some ap- 
proximations in the interpretation of the temperature response data, t*, 
or else devise a means by which the actual response data might be treated 
directly. Some of the difficulties involved with the second alternative 


are: Refer again to equation (9) 
. A 
C7 Ue, pA) = (= o- / — Bur) 
x tis — My. (UH) 
+ Mz (Milvy) € ie POG) LY 


which may be stated as 


N,,. = —by [/- 2204, d+ Ln | / — 9) 


— a 


M. , 


A . 1 CVA) 
+ [Ma elon Ko) C F (/- 9) aV 


and the expression for the inlet fluid temperature, 

Gn) = e776, KF) , 
where @, is the equal to the heater time constant which was derived in 
Appendix C. 

One approach to using equation (9) would be to measure tt at 


" "1 Vv . 
some value of "free time /* > 0 such that gC Ar ) = O for which Neo 


would be, 


NV, —— en Cle & } 
: pt wwala 
+ bn )+ [422 (Mj CVpy) © as “O-¢V)aV 


The difficulty here is, that in order to evaluate 


“M,: ae (NM *(V-)) go a) 1- WV) ALY 


tt rey 
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Nis must first be known. This would result in a direct curve matching 
technique. A second approach would be to measure tt at some time “A70 
such that a ) 

AY) = = [wae x (M2 Cy je 1-9) 

again N. must first be known. In both the above methods if Je is small 
the transient response of the sensing and recording equipment will still 
present problems. 

Because of the difficulties encountered with the above two 
methods, the first alternative, that of making some approximations in the 
interpretation of the measured response data, tts will be considered. 

Since the response of the thermocouples is faster than that of the .003" 
diameter heaters, and approximately equal to that of the .001" heaters, 

see Appendix C, it might be reasonable to assume that the initial response 
recorded during a particular run is that response due to the heaters, see 
Figure 2p, Therefore, it is possible to extrapolate back to "free time" 
Jt = 0, from the time on the response curve where the transient response 
of the heaters has died out and treat the intersection of the extrapolated 
curve and "free time," Lt = O as an actual step change in temperature. The 


above technique is referred to as the "zero intercept" technique. 
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ae, Experimental Technique. 
The existing apparatus at the NPS facility, see Figures 3 and 4, 
has been designed to conform to the idealizations required by Howard [5] 
for the use of his conduction parameter in the maximum slope technique. 
Howard's conduction parameter, ae » is defined as 
i eet das 
where: ne a 
d, = thermal conductivity of the solid (Btu/hr sq ft deg F/ft) 
A, = matrix solid cross-sectional area available for thermal 
conduction (sq ft) 
Howard's idealizations are as follows: 
(1) The fluid flow in the matrix is both steady and uniform 
in velocity and temperature at any cross section 
(2) The matrix thermal conductivity is finite in the direction 
parallel to fluid flow and infinite in the direction normal to flow 
(3) The matrix thermal conductivity is large in comparison to 
that of the contained fluid 
(4) The thermal properties of the fluid and matrix are constant 
and uniform 
(5) The convective heat transfer coefficient is some suitable 
average and remains constant 
(6) A step change in the temperature of the inlet fluid is im- 
posed at real time equal to zero. 
These idealizations result in making the fluid flow one dimensional. 
They in no way conflict with the restrictions imposed by Kohlmayr [11]. 
The requirement of uniform velocity and temperature profiles (1) 
is met by a specially designed entrance nozzle, flow straightening 


screens, and an even distribution of heater wires across the channel. 
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Piersall [15], using the equipment, verified that the velocity and tempera- 
ture profiles were, indeed, uniform within reasonable limits. The im- 
portance of uniform velocity and temperature profiles was demonstrated 

by Wheeler [17]. 

The small temperature change across the matrix of about 20°F gives 
constant thermal properties of the fluid and matrix and a constant heat 
transfer coefficient as required by idealizations (5) and (6). 

The restriction that the temperature of the inlet fluid be subjected 
to a step change will be discussed later. 

The data required for the computation of the various heat transfer, 
fluid flow coefficients, and dimensionless parameters is as follows: 


tee, = atmospheric pressure (mm Hg) 


es = orifice static pressure (inches H,0) 
AR = pressure drop across orifice (inches H,0) 
L\ R = pressure drop across matrix (inches H,0) 


Ps = static inlet pressure at entrance to heaters 
(inches H,0) 


oe = temperature of fluid at orifice (millivolts) 


Che: = orifice diameter 
eG = ratio of orifice diameter to pipe diameter 
Ges = chart speed (sec/inch) 

C3-7, = downstream response 

C3 me? = downstream - upstream response 


Pressures are measured with either a draft gage or a water manometer, 
depending upon the orifice-flow rate combination, except for atmospheric 
pressure which is measured by a mercury barometer. 

The temperature responses recorded are t,-t,> the difference between 
the inlet fluid temperature and the fluid temperature at the matrix exit 
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(used in the maximum slope and zero intercept techniques) and t,-ty» 

the difference between the fluid exit temperature and the matrix inlet 
fluid temperature (used in the centroid technique). For a more complete 
description of equipment, the reader is referred to Appendix A. See 
Figure 5 for the position of temperature and pressure measurements. 

A test run is accomplished by predetermining the necessary pressure 
drop across the orifice to achieve a desired flow rate and to determine 
the number of heaters necessary to achieve a 20°F temperature rise. Air 
is drawn through the apparatus and is controlled at the entrance to the 
turbocompressor. When the desired flow rate is achieved the heaters are 
then energized and the heated air and test core are allowed to reach a 
steady state temperature, at which time pressure measurements are re- 
corded. Power to the heaters is then secured and recordings of t,-t, 


and t.-t both as functions of time are recorded on separate channels 


7 2 
of the Brush recorder (see Figures 2A and B). The temperature at the 
orifice (t) is measured before and after a run to insure the same ambient 
air temperature for the run. 

After completion of the desired runs, the values needed to compute 
slope are taken from the recorded traces of t,-t)- The maximum value of 
slope is obtained visually with the aid of a straight edge. This informa- 
tion is included on the data sheet. The data sheet layout conforms to 
the data input section of a digital computer program [18], which reduces 
the data used in the maximum slope technique and calculates the parameter 


O< used in the centroid technique. For complete details on data re- 


duction the reader is referred to Appendix B at this time. 
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4. Description of Test Matrices. 

The two cores used in this experiment were both cores which had 
been tested previously and for which there was a good deal of data avail- 
able. 

One matrix, the Cercor T20-38, previously tested by Howard [7], is 
a ceramic type-core having a low thermal conductivity, and was used to 
eliminate any errors that might be introduced by longitudinal conduction 
which might have invalidated results using the centroid and zero inter- 
cept techniques. The other core tested, which is a stainless steel plate- 
fin type matrix, Solar 4, had been tested several times [15], [16], [18] 
and was used so that comparisons could be made when longitudinal conduc- 
tion was a factor. 

Further information on core geometries and properties is shown in 


Figures 6 and 7. 
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a. Presentation of Results. 

For each matrix tested, the heat transfer and flow friction character- 
istics have been computed. The computed results are shown tabulated in 
Tables I-II. The values of |» have been plotted for the three methods 
used in this experiment. Figures 8A and 9 compare values of Neu based on 
the maximum slope and centroidal technique for cores T20-38 and Solar 4. 
Figures 8B and 10 compare values of Neu based on the maximum slope and zero 
intercept technique for the same two cores. Figure 11 is a comparison of 
— values computed by the maximum slope technique for the two heater 
systems. Core Solar 4 was the only core tested with the new .001 inch 
diameter wire heaters. 

The heater response curves are shown in Figure 12. See Appendix C 


for complete details used in determining these curves. 
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6, Discussion of Results 

In order to insure that the technique used by the author in this 
experiment was correct, the results of the single-blow test of this experi- 
ment were compared with the results of previous tests of the same cores. 

Considerable difficulty was experienced in the use of the new (.001 
inch diameter) heater system. The heaters burned out frequently, thereby 
limiting the range of Reynolds Numbers over which the core was tested and 
also the number of cores actually tested. The exact reason for the burn- 
out is not known. It was initially thought to be due to surges in cur- 
rent caused by the use of a resistor bank in the first experiments. How- 
ever, when the resistor bank was removed and voltage to the heaters control- 
led by the motor generator set rheostat, the system failed again. In all, 
there were three different failures, all occurring at different flow rates 
and with a different number of heaters in use at the time of failure. 

An examination of Figure 11 shows that in the high Reynolds number 
ranges there is little difference in the value of Neu based on the maximum 
slope data for the two different heater systems. Furthermore, this is in 
a range where the errors in Neu due to errors in the maximum slope techni- 
que are greatest. Without actually applying Kohlmayr's equations in the 
“extension of the maximum slope" technique [17], it is impossible to get 
a quantitative value for the error in Neu due to deviations from the step 
change. However, by the use of his curves based on the "deviation from 
step" (previously defined as I(g) and equal to the area under the inlet- 
temperature response curve) which are presented here as Figure 13, it can 
be seen that for I<~.100 the error in slope, Am which is the difference 
in the maximum slope due to a step change in the inlet fluid minus the 


maximum slope due to a non-step change in the inlet fluid, is less than 
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Pun It must be pointed out that Figure 13, is for an upstream tempera- 
ture change which 1s quadratic, i.e., ECL ) = (1- a). The curves of 
Am versus Be, for the NPS facility would be slightly different from 
those presented in Figure 13, for the g( iL) 's associated with the NPS 
facility's heater: systems are exponential. However, for low values of 
I(I1< .100) the error in maximum slope Am would be small regardless of 
the exact form of the inlet temperature response. Since the maximum I(g) 
resulting from either heater system is less than .03 it would not be 
expected for one to detect any difference of slope and consequently Nie 
resulting from the use of either system. 

Again referring to Kohlmayr's curves, Figure 13, deviations from 
step changes in temperatures would cause large errors at lower flow rates, 
te. Neyo: However, since the I(g) of the NPS facility decreases 
with decreasing flow rates (see Appendix C) it could be assumed that the 
errors in maximum slope due to I(g) would be small. Furthermore, since 
the longitudinal conduction is greater at low flow rates its' effects 
cannot be neglected as they have in Kohlmayr's assumptions. 

In the investigation of the centroid and the zero intercept techni- 
ques the range of flow rates of both cores tested was the same and ranged 
from mM ~ 250 to mM + 950 lbm/hr. The reason for 950 lbm/hr being 
the maximum value tested was that at this flow rate all heaters were in 
use and in order to go to higher flow rates, an increase in voltage would 
have been necéssary to achieve the same 20 deg F. temperature rise. This 
would have changed che heater time constant, which is a function of both 
flow rate and the number of heaters in use. 

Investigation of the curves of Figure 8A and 9 reveals that for both 
cores tested there 1s very close agreement between the results predicted 


by the maximum slope technique, shown as the dotted line, and the results 
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of the centroid technique for values of a between 3.5 and 5. In fact 
there is even some overlapping of points in Figure 8A. For values of 

ee below 3.5, the values of Va based on the maximum slope technique 
begin to decrease from the predicted values whereas the values determined 
using the centroid technique followed closely the predicted values for 
the case of Solar 4 and increased slightly from the predicted values for 
core T20-38. 

It can be seen from Figure 8B, which compares the values of eS 
computed by the maximum slope and zero intercept method, for core T20-38, 
that the flow rates used were not high enough for the zero intercept 
method to be used to its best advantage. Only those runs made at flow 
rates greater than 600 lbm/hr had a clearly measurable zero intercept. 

Figure 10 again compared maximum slope Neu 3 with those evaluated 
by the zero intercept technique in this case for core Solar 4. This was 
done for both sets of heaters. It is important to note the close correla- 


tion between the values of ae for the two heater systems. This will be 


covered more fully in Section 7. 
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7. Experimental Uncertainties. 

Various idealizations and boundary conditions have been imposed for 
the mathematical model of the physical experiment. Due to the fact that 
these idealizations and boundary conditions have not been precisely met, 
certain errors have been introduced. These errors are difficult to as- 
sign a numerical value to, and with the exception of the deviation from 
step temperature change previously mentioned and the effects of longi- 
tudinal conduction, will not be discussed. 

The experimental errors associated with the maximum slope technique 
used at the NPS facility have been covered in considerable detail, [2](3] 
[15]; therefore only those errors introduced in this particular experi- 
ment will be discussed. 

In using Kohlmayr's centroidal technique there are several sources 
of possible error. First, Kohlmayr, in using Hausen's mathematical model 
has neglected longitudinal conduction in the solid. Howard [6] has calcu- 
lated the errors in Neu associated with longitudinal conduction in using 
the maximum slope technique, but no such information is available for use 
with the centroidal technique, for the analytical solution to the single- 
blow problem, in this case, does not include the effects of longitudinal 
conduction. Therefore, one can only get a qualitative idea as to the ef- 
fects of these errors. 

From Howard's curves of Me, versus maximum slope, Figures 14 and 15, 
for a given value of Neu the slope of a temperature response curve is 
different for different values of the conduction parameter x » There- 
fore, the response curves and consequently the centroid of the response 
curves would be different for different values of A. The difference in 
maximum slope due to A decreases as Vig decreases. For example at ao. 


30 a difference in A of approximately .1l causes a 90% difference in 
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maximum slope, whereas for New? the maximum error in maximum slope 
caused by a_ .1 difference in A is approximately 9%. It would seem 
reasonable to assume that provided Kohlmayr's results were applied only 
in the low Ne range Wo» and longitudinal conduction was small (i.e., 
A. <.1), conduction errors could be neglected. For this experiment 

the maximum value of A for either core at values of Neu <_5 was less 
than .022. 

Another source of error is in the determination of the centroid of 
the area under the response curve. Two methods were used for this experi- 
ment, primarily as checks on each other. In both methods the critical 
point is the determination of where the cut-off is for a particular 
curve. The cut-off point is that point on the response curve which is one- 
tenth of the maximum value. This point, defined as Yo occurs at a 
point where the slope of the curve is slowly approaching zero, and an 
error of one millimeter in the ordinate can cause an error of several 
millimeters in the abscissa. For example, in the data taken with core 
Solar 4 when Neu = 1.41, at (t,-t,) = 6mm, O = 13.0 sec., and for (t,-t,) 
= 5, © = 13.75 sec. This was for a typical trace in which (t,-t,)max ~ 
50.5mm. A sketch of the response curve (see mous 16) will help to 
clarify the problem. 

The actual error a caused by an error in detcunimaeyee 
will depend on the value of Neu? A typical example would be for Neu = 
3.5. When computed by the manual technique (see Appendix B for details), 
it was found that a + 1 mm difference in ordinate (which is 7% 2% and 
would be a maximum) caused an error of Y% 2% in enn’ It should be 
pointed out that these values are approximate, for there is some uncertain- 


ty associated with finding the intersection of the three centers of mass 
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lines, Shown in Figure 16 and discussed in Appendix B. This uncertainty 
usually was less than + 2%. The same type of error associated with find- 
ing the cutoff point is also found in the computer technique of determining 
the centroid. In this method the computer solution required data points 
from the response curve at fixed intervals, which for this experiment were 
5mm intervals. Therefore, if the value of AY occurred at a point mid- 
way in any one 5mm increment, there could be as much as a 2% mm error in 
the cut-off. For an average 130 to 140 mm trace this is less than 2%. 

For comparative purposes, data runs were made using both techniques and 
the maximum deviation in centroid between the two methods was found to be 
47%, but the majority of deviations was less than 1%. 

Another difficulty discovered in analyzing the temperature trace was 
caused by the sensitivity of the thermocouples. During a particular run 
a change in the steady flow caused by a sudden draft in the laboratory 
caused a temperature deviation on the trace. This required some visual 
smoothing on the part of the author to extrapolate the actual response. 
This might cause errors when using an automated data reduction process, 
unless it included a good smoothing technique. 

An additional problem was encountered in the process of data reduc- 
tion. Ome of the assumptions made in using the centroid technique is that 
at time zero the temperature both upstream and downstream from the matrix 
is constant. It was found that this was not the case for this experiment, 
the reason being that there was some heating of the thermocouples located 
adjacent to the heaters due to radiation from the heaters. This meant 
that the zero reference point of the trace of (t,-t.) established before 
and after the run (heaters off) was not the same as that just prior to 


time zero when the heaters were deenergized (see Fig. 16). Therefore, 


5 


where the response curve should have resembled curve a of Fig. 16, it 
took the form of curve b. To compensate for this, point c was taken as 
time zero and curve a used to determine the centroid. It should be pointed 
out here that these figures are greatly exaggerated here for clarity and 
the error caused by this difference is felt to be negligible. Also, this 
difference was most evident at low flow rates, for the radiation effect 
was caused primarily by the heaters nearest the thermocouples. As the flow 
rates increased and more heaters were energized the percentage of heating 
due .to radiation became less. 

So far no mention has been made of the error in Bl seg caused by an 
error in Momwrp° Kohlmayr [11] has conducted a linear error analysis 
based on the approximation 

A New ~ tC Mu Dkeenro 
Lhhecenro 

and introduced a relative error amplification factor 


Re wy oh Meu. /Leenrp 


LAcenrp New 


ANtw ~ K AMecento 
Naw JE CENTD 


W sae 
ith empirical results he has tabulated values of K versus Jame and 
K versus N.. 
tu 
° 0 . 

As an example, using Kohlmayr's curves wher 7620 is measured too 
high by 2% and when K = -5.0, then Neu as obtained by this method would 
be 10% too low. 

Associated with the zero intercept technique are two distinct types 


of error. Ome results from the initial assumptions on which this technique 
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is based. These assumptions were discussed in Section 2. If errors are 
present in the results, they may be due to errors in the original assump- 
tions. It is not possible to get, at this time, a quantitative value for 
this type of error or in fact to even verify that the technique used here 
is valid. The only factor which would substantiate the technique's vali- 
dity is that for both heater systems used in the experiment the results 
in terms of Baas were nearly identical. The second type of error is that 
associated with the physical data reduction, primarily the correct extra- 
polation of the response curve back to time zero, see Figure 17. This 


may be seen in the following development: 


- — met) 


Mm, + AN, =— &red-ert At?) 


Me, FON, - 
e , on = = jes Tet, 
e Mer jae 
/—e 


or 


= AN, _ Ya & Ace 
ies 
ej 


wo that if 1-t¢ = .2 corresponding to a value of Neu = 1.609, and 


x 
AT, 
ae «Ao. or 2% 
[- Cy : 
then 
aM, = +. OZO02 or 1.25% 
if AG = ae QZ. or -2% 


/-e 


a 


then 


nko —_ OES Or = lin 
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Gr. Conclusions. 

1. The present design of the new .001 inch diameter heaters was 
not found to be satisfactory due to their frequent and sometimes un- 
explained failures. No improvement of the results of the maximum slope 
and zero intercept tests was obtained using the smaller diameter wire 
heaters. No results were obtained with the new heaters for evaluation 
with the centroid technique, but it should be pointed out that for the 
centroid technique it is not necessary to have a step change in the in- 
let temperature so long as the actual change is known. 

2. It would be extremely difficult to try and compare the three 
techniques used in this experiment on the basis of test results ac- 
curacy. The best that one can do at present is to compare the various 
experimental results with the predicted values based on the maximum 
slope technique evaluated at high values of New (Neu > 5.0). In this 
respect the maximum slope technique is unreliable for Nu 3.5. The 
centroid technique appears to have good results for Neu 5 to as low 
as ./75 which was the minimum value tested in this experiment. For both 
the cores tested using the centroid technique the results were either the 
same as those predicted by the maximum slope technique or somewhat higher. 
The zero intercept method appears to be impractical for values of Neust 2° 
due to the larger differences associated with the logarithm of very small 
numbers. This might be better understood by looking at the governing 


equation in the zero intercept technique 


Nin = — Any (1-£2(Ne,0)) 


and noting that as the argument of the logarithmic term gets smaller the 


natural logarithm itself approaches negative infinity. Below the value of 
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Neu 2? the results were higher than predicted but were extremely con- 
sistent, that is, for different runs conducted at the same conditions 
the results were nearly identical. 

From the standpoint of ease of evaluating, the zero intercept techni- 
que is very fast and simple to apply. Furthermore, in the method used in 
this experiment the zero intercept technique appears to be insensitive 
to slight deviations from a step temperature change. On the other hand 
the results of the centroid technique rely entirely upon how well the so 
called "deviation from step" change of the fluid‘'s upstream temperature 
is known. 

The centroid technique is much more time consuming to use and in 
view of the fact that its results are based on a considerable portion of 
the fluid temperature response trace it is more subject to errors caused 
by sudden fluctuations in the ambient conditions, for in the zero inter- 
cept technique only the very first part of the trace is of interest. 
Based upon the above conclusions the following table is recommended as 
a guide in determining which of the aforementioned techniques should be 


used for a given test range. 


ow 





Neu Range Technique 
oe Neu Maximum Slope 
Leo Ney oe? Centroid Technique 


Ney 2-9 Zero Intercept Technique 





SO eS 6 gE ee ee er a 
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q. Recommendations for Further Study. 

It is recommended that an attempt be made to determine both experi- 
mentally and analytically the effects of deviations from experimental 
assumptions such as constant fluid properties and the convective heat 
transfer coefficient on the results using the different techniques 


including the cyclic testing technique. 
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APPENDIX A 
Description of Equipment 
Air Supply. 

The working fluid for the experiment is air drawn through the test 
equipment by a 30 HP, multi-stage Spencer Turbo-Compressor, which is 
rated at 550 cfm operating on a 220 V a.c. power supply see Figure 3. 
Flow Measuring System. 

Flow measurement was accomplished with an ASME standard orifice 
section. Pressure taps were located d and d/2 diameters upstream and 
downstream respectively from the orifice. Thin concentric orifices with 
throat diameters of .775, 1.232, 1.540, and 2.310 inches were used. 
Heater System. 

Two different heater systems were used. One system utilized a .0031 
inch diameter nichrome wire as the heating element, and the other .001 
inch nichrome wire. Both systems were designed to give the same nominal 
resistance for a given number of heater switches. 

The .0031 inch system, (Figure 18A), consisted of 14 separate bakelite 
frames, each wound with two parallel-connected heater elements. Each 
pair of heaters was controlled by an individual switch, which in turn was 
wired in parallel with the other heater switches. All were controlled 
by one master switch. 

The .0010 inch system (Figure 18B) consisted of the same number of 
bakelite frames, with each frame containing six heater elements connected 
in parallel. The number of heaters and voltage could be varied to achieve 
approximately a 20°F. temperature change for any given flow rate. 

Matrix Holder and Test Section. 


The matrix holder and test section (Figure 4) were both constructed 
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of polyethylene plastic. The matrix holder is a drawer which slides into 
the test section. A removable frame on the side of the holder allows for 
the insertion of the matrix. The flow channel is 3-1/16 inches by 3-1/16 
inches, and can hold matrices up to 3 inches long. The matrix is held 
snugly in place by styrofoam insulation. On the downstream side of the 


holder is a plate containing the thermocouples used to measure t The 


3° 
test section into which the holder slides contains the heaters, pressure 
taps, and the thermocouple set downstream from the heaters, which measures 
=" 
Inlet Cone and Flow Straightener 
This section was designed and tested by Piersall [15], and it pro- 
vided a uniform velocity profile to the air entering the matrix. 
Pressure Measuring System 
Pressure taps are located in the test section upstream and down- 

stream from the matrix holder, (see Figure 5). Two other taps are located 
at the orifice section. Each pressure tap is connected by flexible tub- 
ing to its corresponding manometer and/or draft gage. The following in- 
struments were used interchangeably, depending upon the flow rate: 

1. Ellison Draft Gage Company, 0-3 in. inclined gage 

2. Ellison Draft Gage Company, 10 in. manometer 

3. Ellison Draft Gage Company, 20 in. manometer 

4. Merriman Instrument Company, 120 in. manometer 

5. Precision Thermometer and Instrument Company, mercury 

barometer 

Temperature Measuring System 


Temperatures in the system are measured at four locations; (See 


Figure 5)inlet to the system (t,); downstream from the heaters (t.); 
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downstream from the matrix (t ), and at the orifice (t,) « 
3 


Temperature t, is measured by two different sets of 30 gage iron- 


1 
constantan thermocouples, which were made by Traister [16]. Each set 
consists of 5 thermocouples connected in series. Each thermocouple is 
individually wrapped in teflon tape to prevent shorting. All ten are 
contained in an open-faced aluminum tube mounted in a frame at the exit 
of the inlet cone. The aluminum tube shielded against radiation from the 
heaters. 

Temperature - was measured by a set of five .001 inch diameter 
iron-constantan thermocouples connected in series. The output of these 
was bucked against one of the sets measuring ty so that the output of the 
two sets measured (t,-t,). 


tu 


Temperature t. was measured in the same manner as t, and the output 


bucked against the other set measuring t, so that the output measured 


1 
(t,-t,). The outputs (t,-t,) and (t,-t,) could be recorded separately or 
bucked against each other to give (t,-t)- The desired outputs were then 
led into an Astrodata Model 886 Wideband Differential D. C. Amplifier, 
where they were amplified 100:1. From there the signal was led into a 
Brush Mark 280 Strip Chart Recorder. 


Temperature t,. was measured by a 30 gage copper-constantan thermo- 


0 
couple referenced to an ice junction. The output was read on a Leeds & 
Northrup Millivolt Potentiometer. 
Heater Power. 

The power to the heaters was supplied from the 250 V. D. C. source 


in the laboratory. The actual voltage was controlled both at the supply 


panel and at the apparatus by means of a 50 SL wound resistor. 
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APPENDIX B 
Data Reduction Relationships 
(1) The following is a summary of the data reduction relationships 
used in calculating ws by the maximum slope technique. 
Geometry. 
Three geometric parameters may be used to define compact heat trans- 
fer surfaces. This allows for the comparison of different matrices. 


1. Hydraulic Diameter 


4x free flow area 
Dy , 40), ~ heat transfer area (B-1) 


2. Porosity 


. aree flee auea = re (B-2) 
frontal area Arey 


3. Area Compactness 


B* heat transfer area = A (B-3) 
\ matrix volume A. L 
fr 
Dividing (B-2) by (B-3) 
as P/B (B-4) 
Mass Rate of Fluid Flow 


The mass flow rate, nm, is calculated from ASME Power Test Code [1] 


as modified by Murdock {14] by the following equation: 
e 2 aoa 
Pree ea, 2 A v (B-5) 


where 


flow coefficient including velocity of 


Keroro 
V/ Bt approach 


orofice coefficient of discharge [14] 


oO 
MN 


ratio of orifice diameter to pipe diameter 


D 
Ml 


orifice diameter in inches 


A. 
i] 
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F = thermal expansion factor 


oe = P/RT the specific weight of fluid flowing, assuming a 
perfect gas (1bf/ft3) 


Y = expansion factor 
P = absolute static pressure at orifice (lbf/sq ft) 
R = gas constant for air: 53.35 (ft-lbf/lbm°R) 
T = absolute temperature at orifice (deg R) 

VAs = pressure drop across the orifice in inches H,O 


2 
Substituting the expressions for K and @ in equation (B-5) yields: 


m= 359_f€ BAY [APP 


= aaa (B-6) 
_/ 1-3" RT 


From [1] Fig. 40A 


| — (0.4 + 0. 35S") K 


Y 


k 


1.4 for air, ratio Lp [ko 


Also, from [1], Fig. 38 


F = 1.0 
a 
Z 
= 26) S * ; 
ig Po /13 6)(0.4912*144) Ibf/ft 
eats local atmospheric pressure in inches Hg 
PS = static pressure upstream of the orifice plate in inches Hg 


Making the above substitution in (B-6) with the constants necessary to be 


dimensionally consistent yields: 


m =589,81(Ch//-B" ) a? | 1-(# 


(B-7) 





ea 
at. OSM ag fe | ane lem- @A 2.0) 
(em B72 |e, 8489.7 
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Reynolds Number 


Reynolds Number is defined as: 


Nee = Du G (B-8) 


Ja 


where G, the mass flow velocity = m/A = m/pA, (B-9) 
and Le is the fluid viscosity. 
Substituting: Z 

Mec = ™ OD, [fe Az, P = 4m/ik Fy, B 


or 


W, = #mL [Ob 


(B-10) 


Maximum Slope 


The maximum slope of the downstream cooling curve (t,-¢)) is a unique 


function of N. (13) and A [6]. 


wf ae 
d(== = BIN, , A) 
L(2/rleu) | ma 
oe» Ne and A have been previously defined, but for convenience are 


restated: 
t= FA o Z N, = AA and A= KG As 
Ws & ~ 7x, Mk L 


the new temperature introduced is t. and is equal to the temperature 
of both the solid and the fluid prior to cooling. 


Therefore: 


oN, ail (mk. | Wt, ) Oo 
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and 


LEN.) = (ML/ Wh) LO 


letting; 


(B-11) 


Cc, = MC (Fluid stream thermal capacity rate) 


and 


C. = Wc, (Matrix heat capacity) 


equation (B-11) becomes, 


L(2/M,,) = (Cc/C,) LE 


Furthermore 
ore t, -té 


whose derivative is: 


ol —" = wie ol (ts -t,) (Bap) 
eres eG 


combining equations (B-11) and (B-12) yields; 


d aan 








ere) | Wee 1 ol (ts-t.) 
CAE / Ne NL, em Le (B-13) 


At=t,-t, 


Ge -t, ) 
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From the above sketch of the downstream temperature response curve: 


a (t; =) 





= Y/X 
I ae ' MAX f 
X/ chart speed =~ ao 
d(t,-t,) ae 4 
c=, *At 
combining with the ratio of, 
matrix capacity a (1/sec) 
flow stream capacity rate Ce 
and equation (B-13) 
=e-¢, 
SS) = & al a Chart S$ Peed. 
(B-14) 
dey, ) Cr At X 
| tu 
L MAX 


This value of maximum slope and ,) are then used to enter Table III 
or Figure 14 or 15 to get the corresponding value of Neus 
(2) Centroidal Technique 

For the centroidal technique the recording trace of (t,-t,) is 
used. Two methods are actually involved in utilizing this trace. One 
method is to copy the trace physically with carbon paper on to a thin 
piece of cardboard. Then the trace copy is cut out, and by means of 
a "plumb-line'' attached to a pivot, which allows the trace to swing free, 
the line passing through the pivot point and the center of mass is determined. 
This is done for several points (at least 3) and the intersection of the 


different lines is the centroid of the area under the response curve. 
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These Lines are shown in Figure 16. Ideally all lines will intersect at 
the same place, however, this is not the case in actuality due to errors 
in manipulation. 

An alternate method is also used to find the centroid of the function 
(t,-t,). This involves the use of a computer program, which when fed 
actual curve data (iGmesiNeSa) computes the centroid. This is done in 
the absence of an apparatus which could convert directly the temperature 
response to information useable by the digital computer, such as a paper 
tape puncher. 

Once Maney 18 determined from the trace of (t,-t.) it is necessary 
to compute I(g) in the manner described in Appendix C. This value of I(g) 
is then used to determine LX, » Recall Equation (14) 

P’cemo = fs + .99 TG) + .O/ 
or 

L/S = ae By, L(9) + .Ol 
Once Le is determined, Kohlmayr's curve of Ny vs- Mo » Figure l, 
may then be entered to give Neue 

(3) Zero Intercept Technique. 

To interpret the physical data, i.e., the recorded trace of 


(t,-t))s which represents t._, im the manner described in Section 3, of 


f° 
this report: 

(i) Extrapolate by means of a French curve the response 
curve (t,-t,) from some position on the trace where the transient response 
of the heaters has died out, back to the vertical line passing through the 
point 4 = 0. This point would be the position where the first change in 
temperature was detected by the thermocouples and is referred to as (t4-t,)o 


or ty, (N © ), Figure 2B. 


tu 


(ii) Record the value of (t,-t,) at the intersection of 
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- a 


the time zero line and the extrapolated response curve, also record the 


value of (t,-t)) at the time just prior to the heaters being deenergerized. 


This temperature is referred to as (ty-t) ax 
(iii) Me is then equal to, 


oe = dad! - (t4-ty Jo /(tz-ty) 


which was developed in Section 3, as, 


‘> Lin -t Ne, ,0)) 


(L 


Do 


APPENDIX C 

To determine the transient response of the NPS facility for differ- 
ent heater sets, runs were conducted in which the temperature difference 
(t,-t,) was recorded as a function of time. This temperature difference, 
which was previously called “upstream temperature" is the difference in 
temperature of the air entering the test rig and the air leaving the 
heaters. 

Analytical determination of the heater time constant is straight- 
forward and presents no difficulty based strictly on the analysis of a 
cylinder in cross flow. The analytical determination of the thermocouple 
time constant is more difficult and requires some judgment. Since the 
thermocouples are constructed of .001 inch diameter wire, it is not pos- 
sible to describe precisely and mathematically the geometry of their 
junction. In the construction of the thermocouples the two different 
wires used crossed, at an angle approaching 180 degrees, and arc-welded 
at the junction. Therefore, the junction was neither spherical nor 
cylindrical (which is the case when the wires are butt-welded). For this 
reason the model of the junction was arbitrarily chosen to be spherical 
and calculations were based on a spherical junction .002 inches in dia- ) 
meter. On this basis it was determined that the time constant of the .003 
inch heaters was nearly three times greater than the time constant of the 
thermocouples, while the time constant of the .001 inch heaters was ap- 
proximately two-thirds that of the thermocouples. 

Since the response of the heater-thermocouple circuit would be the 
sum of at least two exponential terms, which would be difficult to 
separate experimentally, the terms have been treated here as one expon- 


ential term, which is a good approximation of the actual system response 
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if one neglects the initial points on the response curve. From the 


normalized temperature response curves the following empirical relations 


were developed: 


where a = time in seconds and 


9/ G, is the heater time constant 
- eo 


FY 
a Gy = in seconds (C-1) 


(4, =>. 


lo] 4-4) ] 2-@ lege (c-2) 
|. -toran Gy 





log € ao 


¢ = 
fy] at 
neue 
(é, 1 ) iy 
where the value of i =e t.) is the reciprocal of 
1 Hee) 
| (t,-t,)max. 


the siope of the temperature response curve. 
The heater time constant was then plotted against mass flow rate (/M). 


From these curves the following relationships were obtained: 


-a logm + log c where asc are (C-4) 


log GO, 


oy = Lin constants (C-5) 


Since the centroid technique makes use of the deviation from step defined 


a IG) = sh GY) Oe (c-6) 


and since JS = ~O6 and Au = xe 
then gY) = 969) = eC. Ce, eo Hao, (C-7) 


Z@) -F Pe ie, du ae AMG, cll 


A= 1] ze e 


where again W/ = /Y)}* constant (for any given 
s&s 


core) 


Therefore to determine I(g) for a particular run, enter the curve of on 
VS. mM (Figure 12) with appropriate mM - Note that this curve is 
independent of the matrix used in the run. Then multiply oy by the 
appropriate @X< ES obtain I(g) used in centroid technique. Figure 12 
shows the experimental values of G,, VS. my plus the theoretical 


plots of heater and thermocouple transient response for both sets of 


heaters. 
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Figure 10, Core Solar Test Results 
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